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ABSTRACT: The average twin lambing rate of Bayanbulak sheep is 2% to 3%. However, a flock of sheep with a close genetic 
relationship and an average of 2 to 3 lambs per birth has been found recently. To determine the major genes controlling the prolificacy of 
the flock in the present study, the flock was designated A while 100 normal Bayanbulak sheep were randomly selected to comprise the 
control flock B. Ligase detection reaction method was applied to detect and analyze the 10 mutational loci of the 3 candidate prolificacy 
genes including bone morphogenetic protein type I receptors, bone morphogenetic protein 15, and growth differentiation factor 9. The 
10 mutational loci are as follows: FecB locus of the BMPR-IB gene; FecX', FecX'^, FecX^, FecX", FecX^, and FecX* of the BMP15 
gene; and Gl, G8, and FecTT of the GDF9 gene. Two mutations including BMPR-IB/FecB and GDF9/G1 were found in Bayanbulak 
sheep. Independence test results of the two flocks demonstrate that the FecB locus has a significant effect on the lambing number of 
Bayanbulak sheep. However, the mutation frequency of the Gl locus in GDF9 is very low. Independence test results demonstrate that 
the GDF9 locus does not have a significant impact on the lambing performance of Bayanbulak sheep. Among the 10 detected loci, 
BMPR-IB/FecB is the major gene that influences the high lambing rate of Bayanbulak sheep. (Key Words: Bayanbulak Sheep, Fertility, 
BMPR-IB, GDF9, BMPI5) 



INTRODUCTION 

Fertility is an important trait of animals, specifically for 
sheep. Three transforming growth members, BMPR-IB, 
BMP15, and GDF9, are the major genes of high prolificacy. 
A746G mutation in exon 6 of the bone morphogenetic 
protein receptor IB gene has been found which induces the 
249th amino acid to turn from glutamine (Q) to arginine (R) 
in Booroola sheep with high prolificacy. The mutation is 
closely linked to the high prolificacy of Booroola Merino 
ewes, and the A746G mutation of the BMPR-IB gene 
located at chromosome 6 has been defined as FecB 
(Mulsant et al., 2001; Wilson et al., 2001). 

By promoting granulosa cell mitosis, bone morphogenetic 
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protein 15 inhibits the expression of the follicle-stimulating 
hormone receptor in granular cells, and adjusts their 
proUferation and differentiation. Thus, bone morphogenetic 
protein 15 plays a key role in the reproductive regulation of 
mammaUan females (Otsuka, 2000; Otsuka, 2001; Juengel, 
2002; Moore and Shimasaki, 2005). The BMP15 gene of 
sheep is located in the X chromosome. Six mutation loci of 
the BMP 15 gene exert a significant effect on the high 
prohficacy of six sheep breeds, namely, FecX' (Inverdale 
sheep), FecX^ (Belclare and Cambridge sheep), FecX'^ 
(Lacaune sheep), FecX^ (Hanna sheep), FecX^ (Belclare 
and Cambridge sheep), and FecX'^ (Rasa Aragonesa sheep) 
(Galloway, 2000; Hanrahan et al., 2004; Bodin et al., 2007; 
Martinez-royo et al., 2008; Monteagudo et al., 2009). 

Growth differentiation factor 9, a growth factor secreted 
by the oocyte, has an important regulatory effect on the 
growth and differentiation of the ovarian folUcle. Sheep 
GDF9 gene is located on chromosome 5. Nine mutations 
{Gl to G8 and FecTT) on the GDF9 gene are related to 
prohficacy (Hanrahan et al., 2004; Nicol et al., 2009). 

Bayanbulak sheep are a local breed which originated at 
the Bayangol Mongol Autonomous Prefecture in Xinjiang, 
China. In addition to good adaptability to the sub-alpine 
steppe climate with an altitude of 3,000 meters, Bayanbulak 
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sheep possess advantageous attributes such as precocity, 
crude feed tolerance, and so on. However, the breeding rate 
of Bayanbulak sheep is 91.5% to 103%, and their twin 
lambing rate is 2% to 3%. In March 2011, we accidentally 
found a family (20 head) of Bayanbulak sheep including 
one progenitor ewe (702), one progenitor ram (707), six 
daughter ewes and one first-descendant ram of the 
progenitor ewe (702), five granddaughter ewes, and four 
daughter ewes and two granddaughter ewes of the 
progenitor ram (707) (Figure 1). The ewes in this flock 
produce two or three lambs per lambing. Their multi- 
lambing rate is 100%. In the current research, the subjects 
are the normal Bayanbulak sheep flock and the small 
Bayanbulak sheep flock with high prolificacy. The 10 loci 
of the 3 genes related to high prohficacy (A746G 
mutational site on BMPR-IB gene exon 6; FecX', FecX'^, 
FecX^, FecX", FecX^, and FecX'^ mutations on the BMP 15 
gene; and Gl, G8, and Fecrr mutations on the GDF9 gene) 
were detected by Ugase detection reaction method (LDR). 
Single nucleotide polymorphisms of the 10 loci were 
detected and compared between the normal flock and the 
high-prolificacy Bayanbulak sheep flock, respectively. The 
current paper aims to determine the major gene responsible 
for high prohficacy in the Bayanbulak sheep flock. 

MATERIALS AND METHOD 

Experimental material 

A small flock of Bayanbulak sheep with close genetic 
relationship, including 20 sheep in the sheep family with 
high prohficacy per lambing owned by a herdsman in 
Bakeshu Village, Hejing County, Xinjiang Province, is 
marked as Flock A (Figure 1). Up to 100 adult sheep with 



normal fertility were randomly selected to form flock B 
from 300 Bayanbulak sheep ewes in Bazhou Livestock 
Breeding Station in Bohu County, Xinjiang Province. All 
blood samples were collected from the jugular vein of the 
experimental sheep. About 2 ml blood sample was collected 
from each animal using an EDTA anticoagulant. The well- 
mixed blood samples were preserved at -20°C. 

Main reagents 

The main reagents used in the present study were 
AxyPrep-96 whole blood genomic DNA (AXYGEN, USA), 
Taq DNA polymerase (Qiagen Hotstar, Germany), Taq 
DNA hgase enzyme (NEB, UK), and dNTP (Promega, 
USA). 

Methods 

DNA extraction: Genomic DNA was extracted 
according to the instructions in the AxyPrep-96 whole 
blood genomic DNA kit, and 0.8% agarose gel 
electrophoresis was used to detect the integrity of the 
extracted genomic DNA samples. 

Primers and probes: According to the sheep BMPR-IB 
gene sequence (AF: 312016), a pair of primers was 
designed to detect the FecB mutation. For the BMP15 gene 
sequence (NM-001114746), four pairs of primers were 
designed to detect the FecX\ FecX", FecX", FecX^ FecX'^, 
and FecX** mutations. According to the GDF9 gene 
sequence (AF:078545), two pairs of primers were designed 
to detect the Gl, G8, and FecTT mutations. Adjacent 
mutations in the same gene can be detected simultaneously, 
so the same primer sequence was designed for FecX", 
FecX', FecX'^, G8, and FecTT. The above-mentioned 
primers were all designed by Primer Premier 5.0. The LDR 
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Figure 1. The pedigree chart of the high-prolificacy flock of Bayanbulak sheep. The gender, age, lambing performance, and FecB 
genotype were marked. The rectangles represent the males, whereas the circles represent the females. 
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Table 1. Information of the PCR primers 



Primers 



Sequence of the primer (5'-3') 



Product 
length (bp) 



Annealing 
temperature (°C) 



FecB up 
FecB low 

Fec}^ up 
FecX^ low 

FecX", FecX', FecX'' up 

FecX", FecX', FecX^ low 

FecX"^ up 

FecX^ low 

Gl up 

Gl low 

GS,Fec7Tup 

G8, FecTTlow 

FecX" up 

FecX'^ low 



GTCGCTATGGGGAAGTTTGGATG 

TGTTGATGAGGCATGAAAACATCTTG 

GCAGCCAAGAGGTAGTGAGG 

ACGAGCCCTCCTCAAGAGA 

GGCAGTATTGCATCGGAAGT 

GATGGCATGATTGGGAGAAT 

TCCAGAACCTTGTCAGTGAGC 

CAGGACTGGGCAATCATACC 

TGAGGCTGAGACTTGGTCCT 

TGTAGAGGTGGCGTCTGTTG 

GAAGCAAATTGCCCAAGACT 

AGGCGTTCTCCTTTCTCCAG 

CTCTGAGACCAAACCGGGTA 

Cl'Gl'CCAAGl'l rrOGGCAAC 



142 
180 
216 
150 
211 
249 

172/189 



53 
53 
53 
53 
53 
53 
53 



probe was designed according to LDR probe design Table 1, and the information of LDR reaction probe is 



principles (Luo et al., 1996). The deletion in mutation locus 
(FecX") requires a specific probe, so only FAM 
modification on the upstream end of primer 5' was carried 
out. The primer information of PCR reaction is shown in 

Table 2. Information of the LDR probes 



shown in Table 2. All the primers and probes were 
synthesized by Shanghai Sangon Biological Engineering 
Co., Ltd. 

PCR-LDR: The total volume of the PCR reaction system 



Probe 



Probe sequences (5'-3') 



Length of 

LDR 
product 
(bp) 



Fei;;6_modify 

FecB_A 

FecBjG 

Ff6'X'^_modify 

FecX^_C 

FecX° _T 

FecX''_modify 

FecX"_C 

FecX"_T 

FecX'_modify 

FecX' _T 

FecX'_A 

Fec"X^_modify 

FecX'' _A 

FecX'' _G 

FecX^_modify 

FecX'^ _G 

FecX" _T 

Gi _modify 

G1_A 

G1_G 

GS-new_modify 

GS-new_C 

GS-new_T 

Fe(;IT_modify 

FecTT_A 

FecTT C 



YTTTTTTrrTrTTTTTTTCATGCCTCArCAACACCGTCT 




FAM 



TTTT' 

P-AGCCATACCGATGTCCGACCTTTTTTTTTTTTTTTTTT-FAM 

TTTTTTTTTTTTTTTTTTTTTTTTrTTTTTTTTTTTTTCAGGCTCGATGGCCAAAACACG 



77 
79 

82 
84 

87 
89 

92 
94 

97 
99 

102 
104 

107 
109 

77 
79 

117 
119 
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was 20 which included 1 |al DNA (50 ng/|il) template, 2 
^1 lOxBuffer, 0.6 nl Mg^"^ (100 mM), 2 nl dNTP (20 mM), 
0.3 \i\ Taq DNA polymerase (5 U/^l), 4 ^l 5xQ-Solution, 
and 0.4 |il primer mixture (5pM). The rest was ddH20. The 
PCR program was as follows: pre-denaturation at 95 °C for 
5 min, 35 cycles of 94°C for 30 s, aimealing at 53°C for 90 
s, extension at 72°C for 60 s, and a final extension at 72°C 
for 7 min. The PCR product was tested using 3% agarose 
gel electrophoresis. The PCR product was diluted with 
isovolumetric ddH20 and used as template in the LDR 
detection. The total volume of the LDR reaction system was 
10 |il including 1 |il lOxBuffer, 1 |il primer and probe 
mixture (0.5 pM/each), 0.05 |il (40U/^1) ligase (Taq DNA 
ligase enzyme system), 1 |li1 PCR product (100 ng/|Lil), and 
6.95 ]jI deionized water. The LDR program was as follows: 
pre-denaturation at 94°C for 2 min, then 35 cycles of 30 s at 
94°C, and 2 min at 50°C. 

Data analysis and genotyping: In the experiment, 1 ^1 
LDR coimection product was mixed with 1 \x, LABI GS 500 
ROX fluorescent molecular weight standards and 1 |il 
deionized formamide sample. The mixture was heated at 
95 °C for denaturation for 2 min, placed into ice for harsh 
cold condition, and then underwent cataphoresis for 2.5 h in 
5% polyacrylamide gel with 5 mol/L urea. GENESCAN"^^ 
672 was used for data collection and lane line correction, as 
well as for size measurement of the migration section and 
modification of the inner molecular weight standard. 
Genemapper was used for data analysis and genotyping. 

Statistical analysis 

Chi-square compatibility test was used to test the 
Hardy-Weinberg equilibrium of allelic genes in the 
mutation locus of each flock. The Chi-square independence 
test was used to determine significance of the mutation loci 
in each flock on the independence status among genotypes. 

RESULTS 

LDR results and gene distribution 

BMPR-IB gene: The LDR products of three different 

lengths were obtained in detecting BMPR-IB A746G 
indicating the existence of the three genotypes. 

In flock A, homozygous mutant genotype (GG) was 
observed 3 times giving a genotype frequency of 15%, 



mutation heterozygous genotype (GA) 14 times giving a 
genotype frequency of 70%, and wild genotype (AA) 3 
times giving a genotype frequency of 15%. Mutation (G) 
gene frequency is 50% and the wild gene (A) frequency is 
50%. 

In flock B, only the wild genotype was observed, so 
wild gene (A) gene frequency is 100%. Mutations in the A 
flock was in Hardy-Weinberg equilibrium. 

The independence test results of flocks A and B of 
Bayanbulak sheep are shown in Table 3. 

Results show that the FecB polymorphism has a 
significant effect on the lambing performance of 
Bayanbulak sheep {x^ = 98.75, p<0.005). 

FecX\ FecX^, FecX", FecX', FecX^ and FecX" loci of 
BMP15 

The LDR results of the six mutation loci at the second 
exon of the BMP15 of Bayanbulak sheep were analyzed. 
The findings indicate that at FecX', FecX", FecX\ FecX\ 
FecX°, and FecX* the lengths of LDR are 92 bp (the 
representative nucleotide is T), 87 bp (the representative 
nucleotide is C), 102 bp (the representative nucleotide is G), 
99 bp (the representative nucleotide is G), 82 bp (the 
representative nucleotide is C), and 189 bp (lack of 17 bp 
deletion), respectively. At each mutation locus, the LDR 
products have only one length, indicating the absence of 
any polymorphism. 

Gl, G8, and FecTT loci of GDF9 

The LDR results of the Gl mutation locus of the first 
exon, the G8 of the second exon, and the FecTT mutation 
were analyzed. Results show that at the Gl mutation locus, 
LDR products of three lengths exist, indicating the 
existence of the three genotypes. At the G8 locus of GDF9, 
one length of 77 bp was found (the representative 
nucleotide is C), dedicating only homozygote genotype. At 
the FecTT locus of GDF9, an LDR product length of 117 bp 
is gained (the representative nucleotide is A), dedicating no 
polymorphism with only the homozygote genotype AA 
present. 

In flock A, there are two genotypes at the Gl point 
mutations locus, e.g., heterozygous genotype AG was 
observed in one sheep giving a genotype frequency of 5%, 
and the wild genotype (GG) 19, with the genotype 
frequency of 95%. Mutation A gene frequency is 2.5%, and 



Table 3. Genotype of the Bayanbulak sheep for the BMPR-IB gene FecB locus 







Genotype (theoretical value) 




df 






GG 


GA 


AA 


Flock A 


3 (0.5) 


14 (2.33) 


3 (17.17) 


20 






Flock B 


0 (2.5) 


0 (11.67) 


100 (85.83) 


100 






T 


3 


14 


103 


120 


2 


98.75 
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Table 4. Genotype of the Bayanbulak sheep for the GDF9 gene Gl locus 





AA 


Genotype (theoretical value) 
AG 


GG 




df 




Flock A 


0 (0.5) 


1 (1.5) 


19(18) 


20 






Flock B 


3 (2.5) 


8 (7.5) 


89 (90) 


100 








3 


9 


108 


120 


2 


0.865 



the wild gene G frequency is 97.5%. 

In flock B, homozygous mutant genotype AA was 
observed 3 times giving a genotype frequency of 3%; 
mutation heterozygous genotype AG 8 times giving a 
genotype frequency of 8%; and the wild genotype GG 89 
times giving a genotype frequency of 89%. Mutation A 
gene frequency is 7%, and the wild gene G frequency is 
93%. 

Gl locus, mutations in the A flock and B flock reached 
the Hardy- Weinberg equilibrium. 

The independence test results of flocks A and B of 
Bayanbulak sheep are shown in Table 4. 

Results show that Gl polymorphism does not affect the 
lambing performance of Bayanbulak sheep (x^ = 0.865, p> 
0.1) 

DISCUSSION 

Thus far, no reports on the prolificacy gene of 
Bayanbulak sheep have been published. In the present paper, 
10 mutation loci of 3 candidate genes of high fertility were 
studied for the first time. Results show that the BMPR-IB 
gene FecB mutation and GDF9 gene Gl mutation are 
present in Bayanbulak sheep. In the BMPR-IB gene FecB 
mutation, the GG or GA ewes produce two or more lambs 
per lambing. The ewes that produce two or three lambs per 
lambing are all GG type or GA type (Figure 1). In the 
control group B, no FecB mutation genotype was observed. 
Thus, the BMPR-IB gene FecB mutation is the major gene 
that influences the high fertility of Bayanbulak sheep. 

According to the study of Fabre et al. (2006) the 
lambing number and ovulation rate are proportional to the 
copy number of the BMPR-IB gene. One copy of Booroola 
gene from the parents increases by 1.5 ova per lambing, 
thus producing 1 more lamb. For homozygous carriers, the 
increase is 3.0 ova, thus producing 1.5 more lambs (Fabre et 
al., 2006). Existing studies have verified the existence of the 
FecB mutation in some world-renowned species with high 
fertility such as Booroola Merino in Australia (Mulsant, 
2001; Wilson, 2001), the Garole in India (Davis et al., 2002), 
Javanese in Indonesia (Davis et al., 2002), and Small-Tailed 
Han Sheep and Hu Sheep in China (Liu et al., 2003; Davis 
et al., 2006; Chu et al., 2007; Chu et al., 2011). The current 
study detected the FecB mutation in Bayanbulak sheep, 
making it the sixth carrier with BMPR-IB gene FecB 



mutation. 

There are many sheep species with single mutations of 
reproductive capacity genes, yet only four have been 
reported to possess multi-locus mutations of reproductive 
capacity genes, including Belclare/Cambridge (Hanrahan et 
al., 2004) with GDF9 and BMP15 mutations, Lacaune with 
FecL and BMP15 locus mutations (Bodin et al., 2002), 
Small-Tailed Han sheep and Hu sheep with BMPR-IB and 
BMP15 locus mutations (Chu et al., 2007), as well as 
Garole with BMPR-IB and GDF9 locus mutations (PoUey et 
al., 2009). Due to the BMPR-IB and GDF9 mutations, 
Bayanbulak sheep have become the sixth sheep species with 
two prolificacy gene mutations. 

The Gl mutation of GDF9 was discovered by Hanrahan 
et al. (2004) in Belclare and Cambridge. This mutation 
leads to amino acid changes owing to the C— >A mutation at 
260 bp of the GDF9. Thus, the 87th amino acid residue 
turns into histidine from arginine. As these mutations either 
happen at the protease cutting site or at the undressed 
albumen of maturation protein, their changes do not induce 
the functional change of GDF9. PoUey et al. (2009) 
reported that Garole sheep exhibit BMPR-IB and GDF9 Gl 
mutations. However, the relation between the mutation and 
the lambing number has not been revealed yet. This paper 
has also discovered that the GDF9 of Bayanbulak sheep has 
a Gl mutation with a very low mutation rate. In flock A, 
only one sheep was heterozygous (a lamb without lambing 
record); the other 19 sheep are wild type. The relevance 
ratio of the control group (flock B) was higher: three are 
mutation homozygous types and eight are mutation 
heterozygous types. The mutation genotype accounts for 
11%. Flock B is a normal flock of Bayanbulak, so its twin 
lambing rate should be the mean value of the group (2% to 
3%). Thus, the Gl mutation locus of GDF9 is not the major 
gene contributing to the large lambing number per lambing. 
Its relation with the twin-lamb lambing performance in 
Bayanbulak sheep requires further study. 

The Small-Tailed Han sheep in China has high 
reproductive capacity with the FecB mutation, originating 
from the Mongolian sheep. The breed was formed after long 
term selection (Li-Xin and Ding-Guo, 2003). Bayanbulak 
sheep also evolved from MongoUan sheep. From the 
perspective of the breed, detecting FecB mutation from the 
flock with high fertility is reasonable. Thus, FecB mutation 
exists in Mongohan sheep with the genetic basis of high 
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fertility, and breeding a line with high fertility is possible. 

CONCLUSIONS 

With the LDR method, two gene mutation loci with 
associated with fertility, namely, FecB of the BMPR-IB and 
Gl (FecG) of the GDF9, were found in Chinese 
Bayanbulak sheep. No mutation in the BMP15 gene was 
discovered. The mutation genotype of Gl of GDF9 has a 
very low ratio within the twinning flock, and the ratio is 
lower than that of the control group. The independence test 
results of flock A and B indicate that Gl of GDF9 does not 
affect the lambing performance of Bayanbulak sheep. Thus, 
Gl is not the major gene affecting the high reproductive 
capacity of the flock. However, the BMPR-IB FecB locus in 
the Bayanbulak flock with large lambing number per 
lambing presents corresponding gene mutations. They 
produce more lambs per lambing and reach a moderate 
polymorphic level. Independence test results of flock A and 
B show that the FecB polymorphic locus has a significant 
effect on the lambing performance of Bayanbulak sheep. 
Thus, this mutation is the major gene affecting the high 
fertility of the flock. Based on these findings, genetic 
resources and marker-assisted breeding technologies can 
accelerate the breeding of Bayanbulak sheep strains with 
high fertility. 
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